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We show that the factorization formula for non-leptonic B decays to two light flavor non-singlet 
mesons derived by Bauer et al. in the context of soft-collinear effective theory is equivalent to the 
corresponding formula in the QCD factorization approach. The apparent numerical differences in 
the analysis of B — > tt-k data performed by these authors, as compared to previous QCD factorization 
analyses, can largely be attributed to the neglect of known perturbative and power corrections. 



The extent to which hadronic decays of B mesons to 
two light hadrons can be computed from first principles 
in QCD has been the subject of many investigations, fol- 
lowing the statement jjj, |2j that the decay amplitudes 
factorize in the limit of very large B-meson mass. Tak- 
ing the final state to consist of two pions, the amplitude 
can be represented schematically in the form 



A(B -> tttt) = F T * $„. + T * $ B * * $„. . (1) 

In this equation, F B ^ n denotes a physical form factor, 
and $s and are the leading-twist light-cone distribu- 
tion amplitudes of the mesons. The quantities T I,n are 
perturbative hard-scattering kernels involving the two 
scales nib (hard) and yroftA (hard-collinear) , which are 
linked to the other elements of the formula by convolu- 
tion integrals (indicated by an asterisk). 

The authors of 3] suggest a factorization formula (their 
Eq. (24)) similar to 0) in the framework of soft-collinear 
effective theory (SCET) and imply that it is conceptu- 
ally different. They also argue that, even at leading or- 
der in the 1/rrib expansion, there may be an additional 
term on the right-hand side of , corresponding to long- 
distance contributions from cc penguins. They do not 
disprove factorization of charm-penguin loops by provid- 
ing a counter-example to factorization; rather, they state 
that they were not able to demonstrate factorization. 

In this Comment, we explain why the formula given in 
is identical in content to the QCD factorization for- 
mula and why non-factorizable charm-penguin con- 
tributions are of higher order in the 1/mb expansion. We 
also point out that the phenomenological analysis of Q 
neglects important perturbative and power corrections 
which are already known. Once these are included, there 
is little room for significant additional contributions to 
the QCD penguin amplitude. 



EQUIVALENCE OF THE SCET AND QCD 
FACTORIZATION FORMULAE 



We first note that the coefficient function of the 
spectator-scattering term can be represented as a convo- 



lution T u = C n * J of hard and hard-collinear coefficient 
functions. The formula quoted in Q follows from by 
rewriting 



(2) 



with Cj w defined as J * <&s * ^tt, and C u defined by 
the decomposition of T 11 shown above. In addition, in 
the SCET form factor C, B7r rather than the physical 
QCD form factor F B ^ n is used. As discussed in 0, this 
implies another rearrangement of this type. 

In general SCET provides a powerful tool to simplify 
factorization proofs (a task that has not yet been com- 
pleted for the case of B — > 7T7t considered here), but 
the resulting QCD factorization formulae can also be ob- 
tained using traditional factorization methods, as is fre- 
quently done in practice. 

The authors of Q entertain the possibility that the 
hard-collinear scale may be non-perturbative, and hence 
they choose not to factorize ( B7T into J * $b * ^tt as 
indicated above. This is a logical possibility in the QCD 
factorization approach. However, we show below that it 
is not supported by theoretical calculations. 

We disagree with Q on the statement that C, B n <C 
pB^ir - 1S a p re di c tion of QCD factorization, whereas the 
SCET treatment suggests C,f x ~ F B ~^ . The factoriza- 
tion formula states that both terms in Q are of the same 
order in 1/mb power counting and that T 1 and T 11 start 
at 0(a°) and 0(a s (-y/wbA)), respectively. However, it 
does not predict the relative size of the two terms, and 
as we explain below, the numerical value of C, B n depends 
on several hadronic input parameters, which are rather 
uncertain at present. 

The authors of 0] point out that the hard-collinear 
kernel J is universal, so that only a single function * 
appears in J3J, which is the same that appears in the 
factorization of the B — > n form factors. This fact is im- 
portant for phenomenology when one opts to treat the 
hard-collinear scale as non-perturbative, but does not by 
itself represent a conceptual difference between the for- 
mulae given in 0,0] and 0. Furthermore, the useful- 
ness of the universality of J is limited to the approxima- 
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FIG. 1: Charm-penguin loop. 



tion where one neglects radiative corrections to the hard- 
scattering kernels C 11 , since only then does the function 
(f^ reduce to a single number 0, 0- In other words, 
a phenomenological treatment of the hard-collinear scale 
as non-perturbative relies on the approximation that the 
kernels are restricted to their tree-level approximations, 
whereas one of the key features of QCD factorization (as 
opposed to naive factorization) is that one can consis- 
tently include radiative corrections. 



II. CHARM-PENGUIN LOOPS 

In Q it is claimed that there is a possible exception to 
factorization from diagrams with charm-quark loops (see 
Figure^]) . The argument is based on the observation that 
when the gluon virtuality is near the cc threshold, q 2 « 
4m 2 ., the non-relativistic scales m c v and m c v 2 become 
important. Since m c v 2 ~ A numerically, this appears to 
introduce a sensitivity to non-perturbative scales without 
power-suppression in l/m&. The authors of Q suggest 
that these diagrams do not factorize, in the sense that 
the long-distance physics at leading-power in the 1/mt 
expansion cannot be factored into form factors or light- 
cone distribution amplitudes. We emphasize that the 
question of non-factorizable cc effects at leading order 
in the heavy-quark expansion is a different issue than 
the one raised in where it is speculated that power 
corrections to the QCD charm-penguin amplitudes may 
be numerically large. 

Factorization statements, be they derived diagrammat- 
ically or with soft-collinear effective theory, always con- 
cern properties of the amplitude in certain asymptotic 
limits, here an expansion in \jm\,, independent of the 
actual size of the expansion parameter. It is therefore 
important to clearly distinguish the issue of factorization 
at leading order in the 1/mf, expansion from the ques- 
tion of whether there are non-factorizable contributions 
which arc formally power-suppressed, but which never- 
theless may be numerically significant. In the following 
discussion, we focus on the question of factorization in 
the formal heavy-quark limit. 

The most intuitive way of understanding why the 
threshold region does not require special treatment is 
based on quark- hadron duality 0- The integration over 
the gluon virtuality in the range < q 2 < m 2 is weighted 
by the pion distribution amplitude, which is smooth over 
the entire integration region. This provides the neces- 



sary smearing of the loop amplitude, which ensures that 
the result is given by a simple partonic calculation up to 
power corrections, in complete analogy with the stan- 
dard justification of the partonic interpretation of in- 
clusive heavy-meson decays, or cross sections in general. 
The smearing also ensures that one can apply the same 
power-counting arguments that demonstrate factoriza- 
tion of other diagrams to charm-penguin diagrams with 
no need to single out the threshold region. In particu- 
lar, in the non-relativistic situation implied in Q there 
is no need to sum Coulomb ladder diagrams, since these 
do not result in large perturbative corrections after the 
integration over q 2 . 

Even without invoking the duality argument, the fact 
that the charm threshold region comprises only a para- 
metrically small portion of the entire integration im- 
plies a phase-space suppression. This fact has been ne- 
glected in the argument of p|. More precisely, writing 
q 2 = iirij, where x denotes the longitudinal momentum 
fraction of the anti-quark in one of the pions, this region 
is Ax ~ v 2 (m c /mi,) 2 or Aa; ~ Km c /m 2 , whichever is 
larger. In order to study the question of whether long- 
distance cc loop effects are of leading order or not, it 
is necessary to decide how the limit to^ — > oo is to be 
taken. If we define the heavy-quark limit by mj, — > oo 
with m c fixed, one may distinguish several possibilities 
such as m c ~ A, m c v 2 ~ A, or m c v 2 3> A. While 
the physics of the threshold region is very different for 
all these cases, they share the common feature that, for 
the purposes of power counting, the charm quark can 
be considered to be a light quark, and the suppression 
of long-distance cc effects has the same origin as that 
for the corresponding diagrams with light-quark loops, 
which implies Aa; ~ \jm 2 . In addition, since in this re- 
gion x ~ m 2 /m 2 : there is a further suppression due to 
the end-point behavior of the pion distribution amplitude 
(which vanishes linearly as x — > 0). 

If we define the heavy-quark limit as mb. c oo with 
the ratio m c /m& fixed, then there is no power suppres- 
sion due to the phase-space or end-point behavior of the 
distribution amplitudes, but the threshold region is per- 
turbative, up to a small non-perturbative contribution of 
order v 2 ■ v 2 (A/(m c v 2 )) 4 — (A/m b ) 4 @. 

We conclude that charm-penguin diagrams factorize at 
leading power in l/irib. The argument for factorization 
remains valid also for more complicated higher-order pen- 
guin graphs whenever the threshold region is phase-space 
(and end-point) suppressed or the charm quark is heavy 
so that perturbation theory is applicable. 



III. VALIDITY OF PERTURBATION THEORY 
AT THE HARD-COLLINEAR SCALE 



In applying the QCD factorization formula to phe- 
nomenology the authors of ,3j treat the hard-collinear 
scale V m bA as non-perturbative, and hence the quantity 
Cj as an unknown phenomenological function. This is 



3 



justified a posteriori following the result of a phenomeno- 
logical fit (on which we comment below). This line of 
argument ignores the fact that perturbation theory at 
scales of order \fm^K ~ m c has been used successfully 
in many important applications in i?-physics, including 
all determinations of \V u i>\ from inclusive 5-decays, and 
studies of the hadronic decay rate of the r lepton. As 
already mentioned, a serious drawback of treating the 
hard-collinear scale as non-perturbative is that it ren- 
ders the factorization approach unpredictive beyond the 
tree approximation, because only the integral over (j* 
and not its functional form can be extracted from mea- 
surements. 

A systematic way to address the question of the pertur- 
bativity of the hard-collinear scale is to calculate higher- 
order corrections in a s to the jet function J in the prod- 
uct T n = C 11 * J, defined as the Wilson coefficient func- 
tion arising in the matching of certain (type-B) SCETj 
current operators onto four-quark operators of SCETn. 
The next-to-leading order terms have been computed re- 
cently and were found to be small [Bj]. Specifically, for 
the case of light pseudoscalar mesons and an asymptotic 
light-cone distribution amplitude <f> w (x) — 6x(l — x), one 
finds that the convolution integrals over the jet function 
give rise to the series 
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a s {pLi) ( (L 2 ) 



1.31(2,) + 1.00 



(3) 

where pn ~ y/rribA is the hard-collinear scale, As is the 
first inverse moment of the £?-meson distribution ampli- 
tude $b(u/, ;U.j), L = ln(mfccj//if ), and (...) denotes an 
average over $s(w,/Xj) with measure duj/to. While the 
precise form of the B-meson distribution amplitude is 
unknown, the fact that u) ~ A ensures that L cannot 
be large, giving a small coefficient to the next-to-leading 
term. For example, using the results of for the mo- 
ments (L 2 ), (L) the coefficient of a s /^ in © is 2.2 ±0.6 
for /i 2 = 0.5 GeVm-f,. There is thus no evidence that per- 
turbation theory cannot be applied at the hard-collinear 
scale. We also note in this context that the power correc- 
tions from the hard-collinear scale are l/mj suppressed 
(and not 1/y/mb) just as those from the hard scale. 

Since the perturbative corrections to the jet function 
are well behaved, the quantity £j v can be factorized and 
expressed in terms of convolution integrals over light-cone 
distribution amplitudes. The question of the numerical 
value of (f™, and whether it is a small contribution to 
the physical form factor F B ~* Tr ( rests on the properties of 
these amplitudes, as well as on other parameters such as 
the strange-quark mass. At leading order in perturbation 
theory, we obtain for (f 77 the result 
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where the notations of j9j have been used. Taking the 
default values and uncertainties of the input parame- 
ters from this reference, and adding errors in quadrature, 



yields £j 77 = 0.016-0.064, which is small compared with 
typical values F B ^ = 0.24-0.30. (F B ^ = C Bjr + (j* 
when hard matching corrections are neglected.) Taking 
some correlated parameter variations so as to reproduce 
the data on B — > tttt decays, scenario S2 of 9] yields the 
somewhat increased value (f 77 = 0.080. To obtain signif- 
icantly larger results would require a very small value of 
the hadronic parameter Xb- While this is a logical pos- 
sibility, a recent QCD sum-rule calculation of As gives 
a value around 0.45 GeV 0, which is in fact somewhat 
larger than the estimate adopted in 0, • 

These estimates are to be compared to the fit result 
Cj 71 = 0.11±0.03 obtained by the authors of who also 
find that the bulk of the B — > tt form factor comes from 
Cf 71 , giving the very small result F B ^ = 0.17 ± 0.02. 
This picture contradicts the QCD sum rules for heavy-to- 
light form factors, in which £j w must be associated with 
a radiative correction ^(j . The preference of the data for 
a smaller B —> tt form factor together with an increase of 
the hard-spectator scattering contribution to the color- 
suppressed tree amplitude a2 has already been discussed 
in |9j, which however did not arrive at a similarly ex- 
treme conclusion. This discrepancy can be traced to a 
few omissions in the calculation of [3j , each of which has 
a minor effect: the absence of radiative corrections, the 
absence of phases in tree amplitudes, the absence of the 
scalar up-penguin amplitude in T c , the use of asymptotic 
wave functions, and finally, a larger value of \V u b\. When 
these effects are taken into account and combined with 
the most recent experimental data, one finds a signifi- 
cantly smaller value of ( B7r and a larger value of F 5 ^^, 
in qualitative agreement with theoretical expectations. 



IV. THE QCD PENGUIN AMPLITUDE 

We now turn to the discussion of the phenomenologi- 
cal analysis of the B — > tttt data performed in 3] . Our 
principal criticism in addition to what has already been 
described concerns the evaluation of the QCD penguin 
amplitude. It may be written as 



P = a 4 + r^ae + (3 3 



-0.09. 



(5) 



where 04 w —0.023 [a°] — 0.002 [a£] represents the vec- 
tor penguin contribution, r^a^ » —0.038 [a°] — 0.014 [a*] 
the l/m& suppressed scalar penguin contribution, and 
/3 3 w —0.011 a power suppressed and rather uncertain 
penguin annihilation term. (The numbers are based on 
the analysis in . Without errors they should be taken 
only for illustration purposes. In particular, we neglected 
all phases, since they are unimportant for the following 
discussion.) 

In the calculation of 3] the lowest order (a®) and some 
of the a s contributions to 04 (those included in the phe- 
nomenological parameters for hard scattering and charm 
penguins) are taken into account. The term r x ag + (3% is 
dropped, because it is power-suppressed. Now while it is 
true that the factorization properties of power-suppressed 
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contributions in general, and scalar penguin contribu- 
tions in particular, have not yet been investigated to all 
orders in perturbation theory, the large tree-level contri- 
bution to r x ae suggests that if one neglects power cor- 
rections entirely (as done in |3() one is certain to obtain 
a poor approximation. We emphasize that this is unre- 
lated to the charm-quark loops discussed above, which 
appear only in the small a s corrections. By dropping the 
scalar penguin amplitude, the authors of [3j are forced to 
erroneously assign the QCD penguin amplitude almost 
entirely to the charm-quark loops. 

There is considerable phcnomcnological evidence that 
the scalar penguin amplitude is in approximate agree- 
ment with our theoretical expectations. The suppres- 
sion of the pseudoscalar-vector and vector-pseudoscalar 
penguin amplitudes relative to the pseudoscalar-pseudo- 
scalar penguin amplitude 0, as well as the pattern of 
drastically different branching fractions for the decay 
modes B — ► rf'^K^ HH, can be attributed directly to 
the different size and sign of the r x a§ term relative to an 



in the QCD penguin amplitude. We are unaware of any 
other theoretical framework that can explain these facts. 
From such studies of penguin dominated -B-decays we are 
therefore led to the conclusion that there is little room 
for extra contributions to the QCD penguin amplitude. 
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